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Photometric Study of KR Aurigae during the High State in 2001 
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Abstract 

We photometrically observed the VY Scl-type cataclysmic variable KR Aurigae after its final rise from 
the fading episode in 2000-2001. Time-resolved observation revealed that the light curve is dominated 
by persistent short-term variation with time-scales of minutes to tens of minutes. On some nights, quasi- 
periodic variations with periods of 10-15 min were observed. No coherent variation was detected. The 
power spectral density of the variation has a power law component (/~^'^^). The temporal properties of 
the short-term variations in KR Aur present additional support for the possibility that flickering in CVs 
may be better understood as a result of self-organized critical state as in black-hole candidates. The light 
curve lacks "super hump" -type signals, which are relatively frequently seen in VY Scl-type systems and 
which are suggested to arise from tidal instability of the accretion disk induced by changing mass-transfer 
rates. The present observation suggests a borderline of superhump excitation in VY Scl-type stars between 
mass ratios q=0A3 (MV Lyr) and q=0.60 (KR Aur). 

Key words: accretion, accretion disks — stars: individual (KR Aurigae) — stars: novae, cataclysmic 
variables — stars: oscillations 
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1. Introduction 

Cataclysmic variables (CVs) are close binary systems 
consisting of a white dwarf and a red dwarf secondary 
transferring matter via the Roche lobe overflow. The re- 
sultant accretion disk becomes thermally stable in sys- 
tems with high mass-transfer rates (M). Such systems 
are called novalikc (NL) stars, because they lack outbursts 
characteristic to dwarf novae [see Osaki (1996) for a re- 
view]. Among NL stars, there exists a small group which 
shows a temporary reduction or cessation of M from the 
secondary. These systems are called VY Scl-type stars or 
anti-dwarf novae (Warner 1995). 

Although accretion disks in NL stars are thermally 
stable, the disk can be tidally unstable (Osaki 1995; 
Osaki 1996). Numerical simulations have shown that 
this instability (tidal instability: Whitehurst 1988) only 
appears below a certain mass ratio {q = M2/M1): q < 
0.25 — 0.33 depending on calculations (Hirose, Osaki 1990; 
Wood et al. 2000; Whitehurst, King 1991; Molnar, 
Kobulnicky 1992; Murray 1998). In recent years, sev- 
eral systems above this stability limit are known to show 
superhumps (Patterson 1999). Since many of them are 
VY Scl-type stars, there has been a theoretical interpre- 
tation that accretion disks can be tidally unstable upon 
variation of M (Murray et al. 2000) even in intermediate q 
systems. Temporary appearance of superhump signals in 
a recent low state of a VY Scl-type star, BZ Cam (Kato, 
Uemura 2001) may support this interpretation. However, 
observations have not yet fully illustrated the upper q limit 
for superhumpers in VY Scl-type stars. KR Aur (g=0.60) 
is an ideal system to examine such a condition, since this 
system has the smallest q among VY Scl-type stars which 



have not been reported to show (or studied for) super- 
humps. 

2. Past Studies of KR Aurigae 

KR Aur is a variable star originally discovered by M. 
Popova in 1960. Popova originally classified the object as 
a 'unusual nova' (cf. Duerbeck 1987). The object has re- 
ceived much attention since Popova, Vitrichenko (1977) 
claimed, from the detection of a large infall velocity 
of 3200 km s^^ in its H/3 line, that the object may 
be a solitary black hole accreting interstellar matter. 
Popova, Vitrichenko (1978) further reported a long-term 
light curve, which could not fit any existing classifi- 
cations of variable stars, a large radial velocity (-1-130 
km s~^), a steep Balmer decrement, and a spectral en- 
ergy distribution approximated by F^=const. Popova, 
Vitrichenko (1978) considered that these features are 
also consistent with an accreting black hole. However, 
Doroshenko et al. (1977) reported that time-scales of con- 
tinuum variation (a few hours to a few tens of seconds) 
and that UBV colors resemble those of CVs, known as 
"novalike" stars at the time of the reporting. 

Doroshenko, Terebizh (1978) studied short-term light 
variation, and detected ~0.4 mag variations with time- 
scales of several minutes. Liller (1980) studied archival 
plates, and obtained the first long-term light curve of 
this object. Liller (1980) showed that KR Aur some- 
times show excursions to faint states (which was also noted 
by Popova, Vitrichenko (1978)), but no clear periodicity 
was found. Popov (1982) reported an ongoing new fading 
episode. 

The nature of the object was revealed with spec- 



T. Kato, R. Ishioka, and M. Uemura 



[Vol. , 



3 

-4—1 



S 



12 
13 
14 
15 
16 



T — I — I — r 



T — I — I — r 



T — I — I — r 



T — I — I — r 



T — I — I — r 



I ' ' ' ' I 

si' 



■! J- M 

T f ■ 



tl_l I I I I I I I I I I I I I I I l__l I I I I I I I I I I I I I U 



49500 50000 50500 51000 51500 

JD - 2400000 



52000 



52500 



4.5 — 
]_ 



-| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 T 

KR Aur J 



il 



Fig. 1. Long-term light variation of KR Aur between 1994 and 2002 drawn from visual observations reported to VSNET. Large 
dots and small dots represent positive and negative (upper limit: object undetected) observations, respectively. Typical errors of 
visual observations are 0.2-0.3 mag. The open circle represents our snapshot CCD observations. The mean epoch of our time-series 
CCD photometry is shown by an arrow. 

troscopic observations by Shatter (1983) and Hutchings 
et al. (1983), who detected a binary motion from radial 
velocity study. The object was thus confirmed to be a 
CV. Williams (1983) showed that the object shows a spec- 
trum typical of a CV with relatively weak emission lines. 
However, since the characteristics of light variations did 
not fit any known subclass of CVs, long-term light vari- 
ations were subsequently studied in detail mainly from 
photographic survey materials (Gotz 1983; Gotz 1984; 
Gotz 1985; Gotz 1986; Gotz 1987; Gotz 1988; Gotz 1989; 
Gotz 1990; Fujino et al. 1987; Liller, Popova 1984; Popov 
et al. 1994; Antov et al. 1996). All of these observations 
showed rather erratic light curves, and occasional fadings 
with time-scales of months to a few years. From these 
characteristics, KR Aur is now considered to belong to 
VY Scl-type novalike CVs, which occasionally show fad- 
ing episodes (cf. Warner, van Citters 1974; Garnavich, 
Szkody 1988; Greiner 1998; Leach et al. 1999). KR Aur 
was also independently selected as an ultraviolet-excess 
object (KUV 06126+2836: Kondo et al. 1984), for which 
Wegner, Swanson (1990) gave a spectroscopic classifica- 
tion as a CV. 

Popov, Antov (1989) reported rapid variations in KR 
Aur. Biryukov, Borisov (1990) further studied this star 
and reported the presence of 25-min periodicity. The light 
curve presented by Biryukov, Borisov (1990) suggests a 
moderate degree of coherence. 
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3. Observations 

The observations were mainly acquired using an unfil- 
tered ST-7E camera attached to 25-cm/30-cm Schmidt- 
Cassegrain telescopes at Kyoto University. A single snap- 
shot observation during a low state was taken on 1998 
December 21, using an unfiltered ST-7 camera attached 
to a 25-cm Schmidt-Cassegrain telescope. All systems give 
magnitudes close to Re- The images were dark-subtracted, 
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Fig. 2. Overall light variation of KR Aur in 2001 
November-December. 

flat-fielded, and analyzed using the Java'^^-based aper- 
ture photometry package developed by one of the au- 
thors (TK). The differential magnitudes of the variable 
were measured against CSC 1889.700 (Tycho-2 magnitude 
y=10.27, B — V^+0.88), whose constancy during the run 
was confirmed by comparison with fainter check stars in 
the same field. The log of observations is summarized 
in table 1. The total number of useful frames was 3139. 
Baryccntric corrections were applied before the following 
analysis. 

4. Results and Discussions 

4.I. Long-Term Variation 

Figure 1 long-term shows light variation of KR Aur 
between 1994 and 2002 drawn from visual observations 
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Table 1. Log of observations. 



Date 




BJD* (start-end) 


Alt 

N ' 


Mag-'- 


brror^ 


T „+ll 

Inst" 


1998 


December 21 


51169.248- 


-51169.248 


3 


6.408 


0.036 


1 


2001 


November 23 


52237.059- 


-52237.244 


428 


3.466 


0.006 


2 


2001 


November 24 


52238.082- 


-52238.111 


69 


3.464 


0.010 


2 


2001 


November 25 


52239.055- 


-52239.172 


270 


3.526 


0.006 


2 


2001 


November 26 


52240.042- 


-52240.177 


252 


3.376 


0.008 


2 


2001 


November 27 


52241.041- 


-52241.105 


142 


3.369 


0.008 


2 


2001 


November 28 


52242.029- 


-52242.152 


259 


3.508 


0.008 


2 


2001 


November 30 


52244.149- 


-52244.352 


355 


3.427 


0.007 


2 


2001 


December 1 


52245.050- 


-52245.126 


102 


3.320 


0.018 


2 


9nni 


December 4 


52248.211- 


-52248.346 
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2001 


December 6 


52250.007- 


-52250.168 


359 


3.465 


0.006 


2 


2001 


December 7 


52251.116- 


-52251.152 


61 


3.285 


0.011 


3 


2001 


December 11 


52255.003- 


-52255.166 


200 


3.455 


0.010 


2 


2001 


December 13 


52257.104- 


-52257.137 


51 


3.460 


0.017 


3 


2001 


December 14 


52258.109- 


-52258.134 


53 


3.614 


0.015 


3 


2001 


December 15 


52259.261- 


-52259.073 


94 


3.308 


0.022 


3 


2001 


December 16 


52260.073- 


-52260.128 


91 


3.229 


0.018 


3 


2001 


December 17 


52261.092- 


-52261.126 


69 


3.272 


0.009 


3 



* BJD-2400000. 
^ Number of frames. 

Averaged magnitude relative to GSC 1889.700. 
^ Standard error of the averaged magnitude. 
II 1: Kyoto (25-cm + ST-7), 2: Kyoto (25-cm + ST-7E). 

3: Kyoto (30-cm + ST-7E). 
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Fig. 4. Representative nightly light curves at different epochs. All light curves show distinct short-term (minutes to tens of minutes) 
quasi-periodic variations. No systematic variation close to the orbital period (0.16280 d) was observed (cf. subsection 4.2). 
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Fig. 3. Period analysis of the November 23— December 17 
data. The upper and lower panels represent the re- 
sults of period analysis with Phase Dispersion Minimization 
(Stellingwerf 1978) and with the CLEAN algorithm (Roberts 
et al. 1987), respectively. The orbital period is marked with 
a tick on the upper panel. No significant coherent periodicity 
was detected at or near the orbital period. 

reported to VSNET Collaboration. ^ Typical errors of 
visual observations are 0.2-0.3 mag, which will not af- 
fect the following analysis. The light curve shows pre- 
dominant low states (durations: several months to two 
years) when the object is fainter than mag 14, and oc- 
casional high states when the object is typically between 
mag 13 and 14. Such a high occurrence of low states is 
relatively rare among VY Scl-type stars (Greiner 1998). 
The long-term behavior of the system most resembles "su- 
perminimum" of MV Lyr (Robinson et al. 1981; Wenzel, 
Fuhrmann 1983; Fuhrmann 1985). At a closer look, the 
system underwent a small dwarf nova- like outburst around 
JD 2450745-2450755. The object then became inactive. 
Before the fully recovery, the object underwent a rather 
complex brightening around JD 2451580-2451630. This 
behavior is quite reminiscent of the fading and recovery 
processes observed in MV Lyr (Pavlenko 1998; Shugarov, 
Pavlenko 1998; Pavlenko, Shugarov 1999). Although such 
behavior is unexpected for variable mass-transfer rates on 
a usual CV (Honeycutt et al. 1994), Leach et al. (1999) 
showed that, in the presence of heating from a hot white 
dwarf, the irradiation on the accretion disk suppresses 
the thermal instability, which can reproduce the observed 
light curve of VY Scl-type systems. From these observa- 
tions, we propose that MV Lyr and KR Aur comprise the 
most "active" subgroup of VY Scl-type stars. 

Figure 2 presents overall light variation of KR Aur 
in 2001 November-December drawn from observations 
in table 1. Although some irregular variation is super- 
imposed, no major "flare" -like brightening (cf. Popova, 
Vitrichenko 1978) was observed. 

^ (http://www.kusastro.kyoto-u.ac.jp/vsnet/ ). 
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Fig. 5. Power spectra of nightly observations. On some 
nights (November 23, 30, December 11), increased power 
around the frequencies 100—150 d~^ (corresponding to the 
periods of 10-15 min) was observed. There was no common 
period between the nightly observations. 



4.2. Superhumps 

Figure 3 presents a period analysis of the November 
23-December 17 data, after subtracting a linear fit to the 
overall light curve. The upper and lower panels repre- 
sent the results of period analysis with Phase Dispersion 
Minimization (Stellingwerf 1978) and with the CLEAN 
algorithm (Roberts et al. 1987), respectively. The orbital 
period is marked with a tick on the upper panel. No sig- 
nificant coherent periodicity was detected at or near the 
orbital period, indicating that no detectable periodic or 
quasi-periodic variations related to orbital modulations or 
superhumps were present (see also subsection 4.3). 

The absence of superhumps makes a clear contrast 
to other VY Scl-type stars with superhumps (TT 
Ari: Udalski (1988); Skillman et al. (1998); Andronov 
et al. (1999); Kraichcva et al. (1999); Stanishcv, 
Kraichcva (2001), MV Lyr: Skillman et al. (1995), V751 
Cyg: Patterson et al. (2001)). This finding suggests that 
the mechanism proposed by Murray et al. (2000) is ineffec- 
tive in the mass ratio of KR Aur. We propose that there is 
a borderline of superhump excitation between mass ratios 
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Table 2. Parameters of VY Scl-Type Stars*. 



Object 


Porb (d) 


PsH (d) 


Ml 


9 


VY Scl 


0.232 




1.22 


0.32 


KR Aur 


0.16280 




0.59 


0.60 


LX Sort 


0.158432 




0.41 


0.88 


BZ Cam 


0.153693 


0.15634 






V794 Aql 


0.1533 








V425 Cas 


0.1496 








"\ C ...1 

VZ Scl 


0. 144622 




1: 


0.7: 


V751 Cyg 


0.144464 


0.1394 






PG 1000+667 


0.144384 








TT Ari 


0.137551 


0.133160 
0.148815 






DW UMat 


0.136607 


0.1330 


0.9 


0.32 


MV Lyr 


0.1329 


0.138 




0.43 


V442 Opht 


0.124330 


0.12090 







* LQ Peg (Kato, Ucmura 1999) is also known as a VY Scl-type star. 
The orbital period has not been reported. 

t SW Sex star (see Hellier 2000 for a recent review; see also 
Thorstensen et al. 1991). A few other SW Sex stars [PX And (Still 
et al. 1995); BH Lyn (Hoard, Szkody 1997)] are claimed to show 
some degree of low/high state transitions (mainly from spectroscopic 
observations). Such a variation may more represent a change in the 
disk state (cf. Groot et al. 2001) rather than VY Scl-type tempo- 
rary reduction or cessation of M from the secondary. Some DQ Her 
stars (intermediate polars) also show temporary fadings (Garnavich, 
Szkody 1988; Hessman 2000), but they are not usually considered 
as VY Scl-type stars. 

References: VY Scl: Martinez-Pals et al. (2000); KR Aur: 
Shatter (1983); LX Ser: Young et al. (1981); BZ Cam: Patterson 
et al. (1996), Kato, Uemura (2001); V794 Aql: Honeycutt, 
Robertson (1998); V425 Cas: see Kato et al. (2001) and the refer- 
ences therein; VZ Scl: Warner, Thackeray (1975), Robinson (1976), 
O'Donoghue et al. (1987), Sherrington et al. (1984); V751 Cyg: 
Patterson et al. (2001); PG 1000-f667: Hillwig et al. (1998); 
TT Ari: Cowley et al. (1975), Thorstensen et al. (1985), 
Udalski (1988), Skillman et al. (1998), Andronov et al. (1999), 
Kraicheva et al. (1999), Stanishev, Kraicheva (2001); DW UMa: 
Shatter et al. (1988), Patterson (1999), Biro (2000); MV Lyr: 
Skillman et al. (1995); V442 Oph: Patterson (1999), Hoard 
et al. (2000), Diaz (2001). 

9=0.43 (MV Lyr) and g=0.60 (KR Aur) in VY Scl-type 
stars. Future determination of mass ratios in longer Porb 
systems (i.e. candidates for systems with higher q than 
MV Lyr) BZ Cam, V751 Cyg and TT Ari is expected to 
provide a more stringent constraint to this limit (see table 
2 for a summary of superhumps and binary parameters of 
VY Scl-type stars). 

J^.3. Short-Term Variations 

Figure 4 shows typical examples of nightly light curves. 
Three longest runs (November 23 and 30, December 
6) were selected as representatives of different epochs 
of the present observation. All light curves show dis- 
tinct short-term (minutes to tens of minutes) quasi- 
periodic variations. No systematic variation close to 
the orbital period (0.16280 d) was observed (cf. sub- 



section 4.2). Aside from the lack of superhumps, these 
variations look remarkably similar to those of quasi- 
periodic oscillations (QPOs) observed in another VY 
Scl-type star, TT Ari (e.g. Mardirossian et al. 1980; 
Jensen et al. 1983; Hollander, van Paradijs 1992; Tremko 
et al. 1996; Andronov et al. 1999). 

Figure 5 shows power spectra of nightly observations. 
On some nights (November 23, 30, December 11), an in- 
creased power around the frequencies 100-150 d~^ (cor- 
responding to the periods of 10-15 niin) was observed. 
There was no common period between the nightly ob- 
servations. This finding confirms the quasi-periodic na- 
ture of the short-term variations. We have not been able 
to confirm the presence of 25-min periodicity claimed by 
Biryukov, Borisov (1990). The periods of the presently 
observed QPOs are close to the periodicities (480-780 s) 
recorded by Singh et al. (1993). We have also confirmed 
night-to-night variation of dominant periods as was orig- 
inally claimed by Singh et al. (1993). 

We also note that the overall profiles of short-term vari- 
ations (figure 4) and the nightly variation of the power 
spectra are also similar to those observed in the peculiar 
symbiotic binary V694 Mon (Michalitsianos et al. 1993; 
Dobrzycka et al. 1996; Ishioka et al. 2001). 

There are two major types of "quasi-periodic" oscilla- 
tions in CVs: dwarf nova oscillations (DNOs), which are 
oscillations observed only in dwarf nova outbursts, hav- 
ing periods of 19-29 s, and have long (several tens to 
^100 wave numbers) coherence times (Robinson 1973; 
Szkody 1976; Hildcbrand et al. 1980), and (in a nar- 
rower sense) QPOs (for a review, see Warner (1995)). 
The present QPOs in KR Aur correspond to the latter 
classification. Several models have been proposed to ac- 
count for QPOs, including vertical or radial oscillations of 
the accretion disk (Kato 1978), reprocessing of the light 
by the orbiting blobs (Patterson 1979), non-radial pul- 
sations of the accretion disk (Papaloizou, Pringle 1978; 
van Horn, Wesemael 1980), radial oscillation of the ac- 
cretion disk (Cox 1981; Blumenthal et al. 1984; Okuda. 
Mineshige 1991; Okuda et al. 1992), excitation of trapped 
oscillations around the discontinuity of physical parame- 
ters (Yamasaki et al. 1995). There has been also a sugges- 
tion that the magnetism of the white dwarf can be respon- 
sible for some kinds of QPOs (Mikolajewski et al. 1990). 
Although present observations are not able to constrain 
the origin of QPOs, a future search for coherent X-ray, 
ultraviolet or optical pulsations, which are a well-known 
signature of a magnetic white dwarf, would be helpful in 
discriminating the possibilities. 

4- 4- Power Spectrum of Flickering 

Figure 6 shows power spectral density (PSD) of the en- 
tire data set. High-frequency white noise has been sub- 
tracted from the PSD. Above the frequency log/=1.5 d^^ 
(corresponding to time-scales shorter than ~45 min). the 
PSD is proportional to f~^'^^. 

These short-term variations (flickering) are one of the 
most characteristic features in CVs. Although the ex- 
istence of flickering in CVs has been well-documented 
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Fig. 6. Power spectral density (PSD) of the entire data set. 
High-frequency white noise has been subtracted from the 
PSD. Above the frequency log/(d~^)=1.5, the PSD is pro- 
portional to /~^'®^. 

since the 1940's (see Bruch 1992 for an extensive his- 
torical review), their physical origin has not been well 
understood. Historically, Warner, Nather (1971) demon- 
strated that flickering almost disappeared during eclipses 
of the eclipsing dwarf nova U Gem. This finding in- 
dicated that the origin of flickering is strongly associ- 
ated with the hot spot (the stream impact point on 
the accretion disk). James (1987) noted the presence of 
power law-type frequency dependence, and proposed that 
a multiple scattering from the flickering source (hot spot) 
is responsible for this frequency dependence. More re- 
cently, Bruch (1991) more extensively studied the proper- 
ties of flickering, and summarized an observational review 
(Bruch 1992). From the analysis of frequency and color 
dependencies, Bruch (1992) and Bruch, Duschl (1993) 
suggested that the inner part of the accretion disk (rather 
than the hot spot) is more responsible for flickering. More 
direct observational evidence for a major contribution 
from the inner accretion disk to flickering has been demon- 
strated through eclipse observations of CVs (HT Cas: 
Welsh, Wood 1995; Bruch 2000, Z Cha: Bruch 1996). It 
is now widely believed that the originally proposed idea 
of stream impact-type flickering (Warner, Nather 1971) 
applies to a only limited sample of CVs, or contributes to 
a small extent to overall flickering activity, and that most 
of CVs have a strong concentration of flickering activity 
toward the inner accretion disk (Bruch 1996). 

In order to reproduce this power-law spectrum, 
Yonehara ct al. (1997) proposed a superposition of "shots" 
in a self-organized critical state (SOC), which was orig- 
inally introduced to explain time variations in black- 
hole candidates (BHCs) (e.g. Mineshige et al. 1994b; 
Mineshige et al. 1994a; Takeuchi et al. 1995; Kawaguchi 
et al. 2000). Taking this analogy into account, the power 
law-type temporal properties of the short-term variations 
in KR Aur present additional support for the possibility 



that flickering in CVs may be better understood as a result 
of SOC as in BHCs. Although the detailed mechanism of 
energy release was not identified at the time of Yonehara 
et al. (1997), WiUiams, Malctesta (2002) recently tried 
to explain flickering by a superposition of flares, resulting 
from an energy release in the photosphere of the accretion 
disk of injected high-energy electrons originating from re- 
connections of magnetic field lines. This type of energy 
release in visual wavelengths would be a promising can- 
didate for explaining flickering in CVs. Further quantita- 
tive comparisons with numerical simulations and observed 
properties in CVs will be a next step toward understand- 
ing flickering in CVs. 

5. Summary 

We photometrically observed the VY Scl-type cata- 
clysmic variable KR Aurigae after its final rise from the 
fading episode in 2000-2001. We show that the long- 
term light curve of KR Aur is densely populated with 
low states, making KR Aur an exceptionally active VY 
Scl-type star. The object showed a complex recovery pro- 
cess from a faint state, which may be understood as the 
result from heating on the accretion disk from a hot white 
dwarf. Time-resolved observation during the high state 
revealed that the light curve is dominated by persistent 
short-term variation with time-scales of minutes to tens of 
minutes. On some nights, quasi-periodic variations with 
periods of 10-15 min were observed. No coherent vari- 
ation was detected. The power spectral density of the 
variation has a power law component {f~^'^'^) at high fre- 
quencies. The temporal properties of the short-term vari- 
ations in KR Aur present additional support for the pos- 
sibility that flickering in CVs may be better understood 
as a result of self-organized critical state as in black-hole 
candidates. Contrary to the large-amplitude short-term 
variations, the light curve lacks "super hump" -type sig- 
nals, which are relatively frequently seen in VY Scl-type 
systems. The present observation suggests a borderline of 
superhump excitation in VY Scl-type stars between mass 
ratios g=0.43 (MV Lyr) and g==0.60 (KR Aur). 
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